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ABSTRACT 

We investigate how a range of physical processes affect the cosmic metal distribution using 
a suite of cosmological, hydrodynamical simulations. Focusing on redshifts z = and 2, 
we study the metallicities and metal mass fractions for stars as well as for the interstellar 
medium (ISM), and several more diffuse gas phases. We vary the cooling rates, star forma- 
tion law, structure of the ISM, properties of galactic winds, feedback from AGN, supernova 
type la time delays, reionization, stellar initial mass function, and cosmology. In all mod- 
els stars and the warm-hot intergalactic medium (WHIM) constitute the dominant repository 
of metals, while for z^2 the ISM is also important. In models with galactic winds, predic- 
tions for the metallicities of the various phases vary at the factor of two level and are broadly 
consistent with observations. The exception is the cold-warm intergalactic medium (IGM), 
whose metallicity varies at the order of magnitude level if the prescription for galactic winds 
is varied, even for a fixed wind energy per unit stellar mass formed, and falls far below the 
observed values if winds are not included. At the other extreme, the metallicity of the intra- 
cluster medium (ICM) is largely insensitive to the presence of galactic winds, indicating that 
its enrichment is regulated by other processes. The mean metallicities of stars (~ Zq), the 
ICM (- 10^1 Zq), and the WHIM (- lO^i Zq) evolve only slowly, while those of the cold 
halo gas and the IGM increase by more than an order of magnitude from z = 5 to 0. Higher 
velocity outflows are more efficient at transporting metals to low densities, but actually predict 
lower metallicities for the cold-warm IGM since the winds shock-heat the gas to high temper- 
atures, thereby increasing the fraction of the metals residing in, but not the metallicity of, the 
WHIM. Besides galactic winds driven by feedback from star formation, the metal distribution 
is most sensitive to the inclusion of metal-line cooling and feedback from AGN. We conclude 
that observations of the metallicity of the low-density IGM have the potential to constrain the 
poorly understood feedback processes that are central to current models of the formation and 
evolution of galaxies. 

Key words: cosmology: theory — galaxies: formation — galaxies: abundances — intergalac- 
tic medium 



1 INTRODUCTION 

The spatial distribution of elements synthesised in stars (henceforth 
'metals') provides an archaeological record of past star formation 
activity and of the various energetic phenomena that stirred and 
mixed these metals. Recent cosmological simulations of galaxy for- 
mation follow the different stellar evolutionary channels through 
which metals are produced, and include some processes that cause 
metals to escape from their parent galaxy. Here we will investigate 
to what extent these simulations produce realistic enrichment pat- 
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terns, which physical processes affect the metal distribution, and 
how robust the predictions are. 

Observations can constrain integrated stellar and neb- 
ular (ISM) abundances of galaxies (for z ^ 1 see e.g. 
Kobulnickv & Zar itskv 1999; Kunth & 6stliij200(]|:lTremonti et al.1 



Dunne et al.li200i f or z 1 see 



Churchill et al. 



200l for 



lErbet al Jl20^ : iMannucci et al j |2009h . Metals in cold 



2004 

z > 2 see 

gas can be observed in absorption against a backgroun d source, see 
the dampe d Lyman-al pha observations at z ~ 3 (e.g. ^Pettini et alj 
J994; Proc haska et aU.2 003) a nd at lower (e g.lMeirhig et al. 20 0^ 
and higher redshift (e.g. lAndo et alj l2007l : IPrice et alj 1200% . If 
the background source is sufficiently bright, then this technique 
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can a l so be a pplied to more di ffuse gas (e.g . Songaila & Cowie 



1996 



200C, 



Cowie & Songaila 199 11: lEUison et all l2000inSchaye et al 



2003i. , 2007.: .Simcoe et ai.l |2004 IScannapieco et al I l2006l : 

Aguirre et al.l2008h . The metallicity of the ICM is inferred from X- 
I '' II 1 I II ) 



ray observations (e.g. Mushotzky et al. 1996; de Plaa et al. 2006': 
Sato et al.l l2007l: [kasm ussen & Po nman 2007: Ma ughan et al. 



20081: Leccardi & Molendi 2008: Snowden et al.ll2008 



These observations use different tools to observe a variety of 
elements in a variety of ionization states, and it is not always ob- 
vious how to convert all these measurements to a common 'metal- 
licity'. Often this is done assuming that the relative abundances of 
elements equal those measured in the Sun. Unfortunately even the 
assumed metallicity of the Sun itself varies between authors. Here 
we assume Zq — 0.0127, and solar abundances from the default 
settings of CLOUDY (version 07.02, last described by Ferland et al. 
1998). 

Convolving metallicity with the fraction of baryonic mass in 
each of the different phases then yields a census of cosmic metals 
(e.g. Fukugita et al. 1998). It must, however, be kept in mind that 
a large fraction of the metals are potentially not accounted for by 
the current data. For example, a large portion of z = baryons are 
thought to reside in the warm-hot intergalactic medium (WHIM), 
which has not yet been convincingly detected. If, as is likely, the 
WHIM is also metal enriched, then it could harbour a significant 
amount of metals. Similarly, the prevalence of hot (T > 10'' K) 
metals at z > 2 is currently poorly constrained. Gas cooling, star 
formation, galaxy interactions, ram pressure stripping, stellar and 
galactic winds all cause gas and hence metals to be cycled between 
the different phases, which makes it complicated to relate the ob- 
served metal distribution to the original source and/or enrichment 
process. 

A theoretical calculation of the census of cosmic met- 
als must take into account their production by nucleosynthe- 
sis, their initial distribution, and the mixing that occurs in later 
stages. Numerical simulations that include such 'chemo-dynamics' 
have become incr easingly sophisticated since the early work by 
lTheisetalJl|l992h . usually concentrating on the evolution of a sin- 
gle g a laxv jSte inmetz & Muller 1995; Raiteri et al. 1996; Berczik 



19991: 'Recchi et al. 2001 : Kawata,2 001: Kawata & Gibson i2003l : 
Kobavashi 2004; Mori & Umemur a 2006; Stins onet alj | 200( 



Governato et al. 2007|; IBekki & Chibai i2009i : iRahimi et al. 20U 
Pasetto et al. 2010), or a cluster of g alaxies (e.g. iLia et al.ll200^ 



Vald arnini 2003; Tornato re et alj20p 4; Sommer-Lar sen et alj2005l : 

"20101) 



Romeo et alji2006 : i Tornatore et al.ll2007i : ,Fabian et al. 



cosmological context. These simulations include interpolation ta- 
bles from stellar evolution calculations for the production of met- 
als, but implement the enrichment processes explicitly, for ex- 
ample by injecting metal-enriched gas near a site of star for- 
mation in a galactic wind. The fluxes of metals between phases 
due to gas cooling and galaxy/gas interactions are computed ex- 
plicitly by these hydro-codes. Early simulations that were also 
used t o look at the met a ls outs ide g alaxies include iCen & Ostriked 
( Il999l) , iMosconi et all ( 1200 ih and iTheuns et al] ( |2002|) . Subse^ 
quent authors implemented more physics while increasing res- 
olution and/or box size in order to minimise numerical ef- 

2005": 



fects (e.g. IScannapieco et al. 



Oppenheinier&Dav 



Cen & OstrikeJ l2006l: iBrook et all i2007,: .Kobavashi et al 



Qppenheimer & Davd200^ : IWiersma et alj2009l.l201ol : IShen et al.1 
l1|2I 



20061: 



20071: 



20091 : iTomatore et al 120091) . 

In this paper we use cosmological hydrodynamical simula- 
tions of the formation of galaxies to attempt to answer the question 
'Wliere are the metalsT , by computing the evolution of the frac- 



tion of metals in stars and variou s gas phases. Our s uite of numeri- 
cal simulations (the OWLS suite: ISchave et alj201(ll) , includes runs 
that differ in terms of their resolution, input physics, and numerical 
implementation of physical processes. We use it to investigate what 
physical processes are most important, how reliable the predictions 
are, and to what exte nt these depend on the sometim es poorly un- 
derstood physics (see Sommer-Larsen & Fvnbdl2008l , for a similar 
invesitgation). In lWiersma et alj j2009l) we introduced the method 
and described some of the numerical issues involved, but analysed 
only a single physical model (the OWLS reference model). Here 
we vary the cooling rates, star formation law, structure of the ISM, 
properties of galactic winds, feedback from AGN, supernova type 
la time delays, reionization, stellar initial mass function, and cos- 
mology. 

This paper is organised as follows. Section |2] introduces the 
simulations used. Those familiar with the OWLS project may skip 
to the results in Section [3] which begins with an overview before 
discussing the physical variations which are the most relevant to 
the metal distribution. This section closes with a summary of all 
the simulations in the OWLS suite in Section [Tsl In Section |4] we 
present our conclusions. 



2 SIMULATIONS 

The OWLS suite jSchave et alj|201oh consists of more than fifty 
large, cosmological, gas-dynamical simulations in periodic boxes, 
performed using th e A''-body Tree-PM, SPH code GADGET III; see 
Tables 1 and 2 in ISchave et all feoid) for a full list of parame- 
ters. GADGET III is an updated version of gadget ii dSpringeil 
|20() ^, to which we ad ded new physics modules for star formation 
( Scha ve & Dalla Vecc hia 2008), feedb ack from supemova e (SNe) 
in the form of galactic winds (Dalla Vecchia & Schave 2008), feed- 
back from accreting black holes ( Booth & Schave 2009) , radia- 
tive cooling and heating in the presence of a n ionizing backgroun d 
dWiersma etalllTooa) . and stellar evolution jWiersma et al.ll2009h . 
In this simulation suite, numerical and poorly known physical pa- 
rameters are varied with respect to a 'reference' model, to assess 
which conclusions are robust, and which processes dominate. The 
cosmological parameters of the reference model, its subgrid and 
their models and their their numerical imp lementation, are d is- 
cussed briefly in the next section. GIMIC jCrain et d] |2009h IS 
a complementary suite of simulations performed with the same 
simulation code, but employs a single set of parameters to investi- 
gate how star formation depends on environment, by simulating re- 
gions picked from a large volume simulation. Here we give a short 
overview of the code, concentrating especially on those processes 
that are directly relevant to metal enrichment. 



2.1 The REFERENCE model 

The REFERENCE model assumes a cosmologically flat, vac- 
uum energy dominated ACDM universe with cosmological pa- 
rameters [Q„,,nb,Q.K,Oii,ns,h] = [0.238,0.0418, 0.762,0.74, 
0.951j0.73], as determined from the WMAP 3-year data and con- 
sistenij with the WMAP 5-year data. The assumed primordial he- 
lium mass fraction is y = Xhc = 0.248. We used CMBFAST (ver- 
sion 4.1: [Seliak & Zaldarriaga.1996) to compute the linear power- 



^ Our value of erg is 8 % lower than the best-fit WMAP 7-year data 
)jarosiketalj201(]|) . 
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spectrum at the starting redshift z = 127. Simulations with given 
box size use identical initial conditions (phases and amplitudes of 
the Gaussian density field), enabling us to investigate in detail the 
effects of the imposed physics on the forming galaxies and the in- 
tergalactic medium. The simulations are performed with a gravi- 
tational softening that is constant in comoving variables down to 
z = 2.91, below which we switch to a softening that is constant in 
proper units. This is done because we expect two-body scattering to 
be less important at late times, when haloes contain more particles. 

• Gas cooling and photoionization Radiative cooli n g and heat- 
ing are implemented as described in 'Wiers ma et al] |200#I . In 
brief, the radiative rates are computed element-by-element, in the 
presence of an imposed ionizing background and the CMB. We 
use the redshift-depen dent ioniz ing background due to galaxies and 
quasars computed by iHaardT ife Madau (2001, hereafter HMOl). 
Contributions to cooling and heating of eleven elements (hydrogen, 
helium, carbon, nitrogen, oxygen, neon, magnesium, silicon, sul- 
phur, calcium, and iron) are tabulated as a function of density, tem- 
perature and redshift, using the publicly available photo-ionization 
package CLOUDY, last described by Ferland et al. ( 1998), assuming 
the gas to be optically thin and in (photo-)ionization equilibrium. 

Hydrogen reionization is implemented by switching on the 
evolving, uniform ionizing background at redshift z — 9. Prior to 
this redshift the cooling rates are computed using the CMB and a 
photo-dissociating background which we obtain by cutting off the 
z — 9 HMOl spectrum above 1 Ryd, which suppresses H2 forma- 
tion and cooling at all redshifts (we do not resolve haloes in which 
Pop. Ill stars that form by H2 cooling). 

• Star formation The interstellar medium of the Milky Way con- 
sists of multiple 'phases': a warm component in which hot super 
nova bubbles envelope and penetrate cold 'clouds'. Current cosmo- 
logical galaxy formation simulations cannot resolve such a multi- 
phase ISM, and in addition not all the physics that governs the in- 
teraction between these phases, and the star formation in them, is 
included. Instead of trying to simulate these physical processes, we 
use the following 'sub-grid' model. ISM gas is assumed to follow 
a pressure-density relation 



P = Po(wH/nH;thr. 



(1) 



where Po /k = 1.08 x 10^ cm"^ K. We use 7off = 4/3 for which 
both the Jeans mass and the ratio of the Jeans length to the SPH ker- 
nel are independent of the density, thus preventing spurious frag- 
mentation due to a lack of numerical resolution. Finally, only gas 
dense enough to be gravo- thermally unstable, wh > nHithros = 
10~^ cm~^, is assumed to be multiphase and star-forming gas 
dSchavdliooi) . It is on this gas that the pressure-density relation 
(Eq.lTJ is imposed. 

Star formation in d isk galaxies is ob served to follow a power-law 
'Keimicutt-Schmidt' jKermicutJ 1 9981) relation between the surface 
density of star formation, E,, and the gas surface density. Eg, 



E. = 1.5 X lO~'*M0yr~^kpc~ 



1 M0PC- 



(2) 



We have implemented this surface density law as a pressure 
law. A star formation rate which depends on ISM pressure guar- 
antees that simulated disk galaxies in which the disk is verti- 
cally in approximate hydrostatic equilibrium follow the observed 



^ We used their Eq. (3) rather than (4) and CLOUDY version 05.07 rather 
than 07.02. 
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Figure 1. Mass weighted metal distribution in temperature-density space 
at z = in simulation REFERENCE_l100n512. The colour scale gives 
the metallicity. The contours indicate the metal mass distribution and aie 
logarithmically spaced by 1 dex. The straight lines indicate the adopted 
division of the gas into: star-forming gas (i.e., nn > O.l cm~^), diffuse 
IGM (p < 10^ (p), T < 10^ K), cold halo gas (p > 10^ (p), T < 
10^ K), WHIM (10^ K < T < 10^ K), and ICM (T > lO'^ K). The 
metals are distributed over a wide range of densities and temperatures. 



law, independent of the value of jctt imposed on their ISM gas 
dSchave & Dalla Vecchiall2008h . 

• Gal actic winds Galact i c wind s are implemented as described 
in Dalla Vecchia & Schavd ( I2OO8I) . Briefly, after a short delay of 
tsN = 3 X 10^ yr, corresponding to the maximum lifetime of 
stars that end their lives as core-collapse SNe, newly-formed star 
particles inject kinetic energy into their surroundings by kicking a 
fraction of their neighbouring gas particles in a random direction, 
as governed by. 



1 2 



(3) 



The mass loading factor, rj, relates the rate at which mass is 
launched in a wind, M^, to the star formation rate, Al,. The prod- 
uct 77 w^, where 77 is the mass loading factor and is the initial 
wind velocity, is proportional to the fraction of the SN energy 
produced per unit mass, esN, that powers the wind. We usually 
characterise the wind implementation by the mass loading factor, 
and wind speed. 

The wind prescription is implemented as follows: each SPH 
neighbour i of a newly-formed star particle j has a probability of 
V^j / '^'^-i" of receiving a kick with a velocity ii„. Here, 
the sum is over the A^ngb = 48 SPH neighbours of particle j. 
The REFERENCE simulation described below has a mass loading 
of = 2 and wind speed of Uw = 600 kms~^ (i.e., if all baryonic 
particles had equal mass, each newly formed star particle would 
kick, on average, two of its neighbours, increasing their velocity 
by Uw ~ 600 kiTis~^). Assuming that each star with initial mass 
in the range 6 — 100 M© injects 10^^ ergs of kinetic energy as it 
undergoes a core collapse SN, these parameters imply that the total 
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wind energy accounts for 40 per cent of the available kinetic en- 
ergy for a Chabrier IMF and a stellar mass range 0.1 — 100 Mq 
(if we consider only stars in the mass range 8 — 100 Mq for type 
II SNe, this works out to be 60 per cent). The value rj — 2 was 
chosen t o roughly repro duce the peak in the cosmic star formation 
rate (Schave et al.ll2O10l) . Note that contrary to t he wid ely-used ki- 
netic feedback recipe of ISpringel & Hernquisll j2003h , the kinetic 
energy is injected locally and the wind parti cles are not decou- 

S iled h ydrodynamically. As discussed by Dall a Vecchia & Schavd 
200i), these differences have important consequences. 
• Chemodynamics We employ the method outlined in 
IWiersma et al. Briefly, a star particle forms with the 

elemental abundance of its parent gas particle. It then represents 
a single stellar population (SSP) with given abundance, and an 
ass umed stel l ar IM F. The reference model uses the IMF proposed 
bv lChabrieJ j2003h , with mass limits of O.IM© and lOOM©. At 
each time step, we compute the timed release of elements and 
energy from three stellar evolution channels: (i) core collapse SNe, 
(ii) type I SNe and (iii) asymptotic giant branch (AGB) stars. The 
stellar evolution prescriptions are based on the Padova models, 
using stellar lifetimes computed in IPortinari et alj 1 1998[) an d the 
yields of low-mass a nd high-mass stars of iMarigol j200lh and 
IPortinari et alj jl998 ^ respectively. The vields of IPortinari et al.l 
1 1998h include ejecta from core collapse SNe (SNII) along with 
their calculations of mass loss from high-mass stars. Using these 
yields gives an element of consistency between the high- and low- 
mass stellar evolution. An e-folding delay time is used to describe 
the SNIa rate. The observed cosmic SNIa rate is approximately 
reproduced in the REFERENCE mo del with a fract i on of 2.5% of 
white dwarfs that become SNIa, see IWiersmaetaI] ( l2009h . We use 
the 'W7' SNIa yields of iThielemann et alj ilOO^ 

Elements produced by nucleo-synthesis are distributed to SPH 
particles neig hbouring the star, w eighted by the SPH kernel, as 
done by e.g., Mosconi et alj j200lh . The simulations track the abun- 
dance of the eleven individual elements that are important for the 
cooling and we use an extra 'met allicity' variable to tr ack the to- 
tal metal mass of each particle (see lWiersma et alj2009l) . The ratio 
Z = Mz/M of metal mass over total mass is the 'particle metal- 
licity' . For gas particles we also compute a 'smoothed metallic- 
ity' as Zsm = pz/p, i.e. the ratio of the metal density, pz, over 
the gas density, p, where both densities are computed using SPH 
interpolation. Stars inherit metal abundances of their parent gas 
particle: we record b oth the particle and smoothed metallicity. In 
IWiersma et al. I j2009l) we argued that the smoothed metallicity is 
more consistent with the SPH formalism than the particle metallic- 
ity. Using smoothed metallicities results in a spreading of metals 
over slightly greater volumes. 



2.2 The Simulation Suite 

Table [T] contains an overview the models from the OWLS 
suite that we consider here. Simulation names contain a string 
'LXXXNYYY' which specifies the co-moving size of the periodic 
box, L=XXXh^^ Mpc, and the number of particles N=YYY'^ of 
(initial) gas and dark matter; most runs discussed here have L=25 
or 100, and N=512'^. Comparing different L and N models allows 
us to investigate the effects of numerical resolution, and missing 
large-scale power. The table also contains a brief description of 
the physics in which a p a rticula r run differs from the REFERENCE 
model, see ISchave et al.l feOlOl) and the next section for more de- 
tails. 



3 RESULTS 



3.1 The reference model 



We begin by describing the metal distribution of the REFER- 
ENCEX100n512 simulation in temperature-density space at red- 
shift z — shown in Fig.[T] 

It is useful to divide the T — p plane into several 'phases', 
since, as we will show, different enrichment processes dominate in 
different phases, and also because the observational constraints for 
different phases are inferred from different types of data. A first di- 
vision distinguishes between star-forming gas (henceforth SF gas) 
on the imposed P ~ p relation (Eq. [T), which we identify with 
the ISM in galaxies, and non-star-forming gas (NSF gas). The SF 
gas can be seen as the thin band of contours on the right of figure 
[T]at p/{p) > lO''. We further divide NSF gas into hoi gas, typi- 
cally found in halos of large groups or clusters (ICM, T > lO'^ K), 
warm-hot intergalactic and circum- galactic gas (WHIM, 10 ' K < 
T < 10^ K), dijfuse gas (diffuse IGM, p < 10^ (p), T < 10^ K), 
and cold halo gas (p > 10'^{p), T < 10^ K). These phases are 
indicated in Fig.[T] The colour coding indicates the mass-weighted 
mean metallicity for a given temperature and density. The metal- 
licity shows a clear gradient with density, with the exception of 
the track of photo-ionised gas in the lower left hand comer, which 
has a very low metallicity. The contours, on the other hand, show 
the distribution of total metal mass. Metals are found at virtually 
all temperatures and densities, from the high-densities in galactic 
disks, to extremely low densities. 

The spatial distribution of metals is shown at z — 4, 2 and 0, 
in Fig.|2| At z = 4 (top left) metals are strongly clustered around 
haloes, with circum-halo metallicities of log{Z/ZQ) ~ —3 to - 
2, and large fractions of volume are enriched to exceedingly low 
levels, or not at all. As time progresses, circum-halo metallicities 
increase to XogiZ j Zq) ~ —2 to -I and the filling factor of metals 
also increases substantially, yet even at z = (bottom right) there 
are still co-moving volumes which are barely enriched. Comparing 
the two resolutions (and box sizes), we note that although the abun- 
dance of circum-halo gas is similar at 2; = 2 in the L100N512 (bot- 
tom left) and L025N5 12 (top right) runs, the filling factor of metals 
is higher in the higher resolution simulation. Indeed, we showed in 
I Wiersma et al.l (^010) that lower density gas (which accounts for 
larger filling factors) was typically enriched by lower-mass galax- 
ies, and higher resolution simulations probe the galaxy mass func- 
tion down to lower masses. 

We continue our analysis by comparing the metallicities and 
metal mass fractions in different phases provided by the reference 
simulation with observations. We then turn our attention to a com- 
parison of the different OWLS models. 

Fig.|3]shows the evolution of the metallicitjlfl metal mass frac- 
tion, and Q,z for various phases in the REFERENCEx100n512 
simulation. The fraction of metals in a given phase is simply the 
ratio of the metal mass in that phase over the total metal mass 
in the simulation. The amount of metals in a phase is shown as 
Jlz = (pz) / Peril, the ratio of the mean metal mass density in a 
given phase over the critical density. As the metal mass fraction 
and Qz are much harder to measure observationally, and are in any 
case inferred from metallicity measurements, we only compare the 



^ Note that for the remainder of this paper, we use 'pailicle metallicities' 
(see section im . although since we are considering mean metallicities of a 
given phase, the results would have been virtually identical if we had used 
smoothed metallicities. 
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Table 1. Simulation Set 



Name 


Box Size (Mpc/h) 


Description 


AGN 


25 


100 


IncoiTDorates AGN model ot Booth & Scliave (2009). 


DRl TMFrONTSFVIfilS 


25 


100 


TV\n-hpnw TA^T-^ for Ti tt ^ ^0 f*vn ' — 1 1 S Vtyi q 


DRT TMFVl (SI 8 


25 


100 


Tnn-hpnvv TMF fnr titt ^ "iO rm~'^- n — 1 fil Vm ci"^ Y ('0^1 — 9nS'^v1 n~^M^vr~^lrr»r~'^ 

1 tJU-ilCaV y liVll iUl /t ^ \J<J f^LLl , i^w XU iO Jvlli o , Z-i * I VJ f — ^ . UOO J. W iVJ.(V) V 1 IVL/L, 


LJUL^LlVll v^Wi"! 1 vJl iVll_j 1 


25 


100 


Tnn-Vipftw TAAF fnr n tt ^ "^O nm n — 1 A '^A'^ 


F\TJT Tl\yTTJl\yTT 1 A 
UljLliVir'iViL 14 


25 


100 


lop-neavy livir' lor tlh s> oU cm , tj — i^.d^o, Zj*i^uj — z.uoo x iu iviQyr Kpc 


lLKJiJ ipu 


25 


100 


Isothermal ec[uatiori of state 


FOS 1 


25 




Emjation of state p oc p'^* , '"y* ^ 5/3 


TA/TF<i AT P 


25 


100 


Salpeter IIVIF, SF law rescaled 


MJi T 
iVllLiL 


25 


100 








100 


AGB & SNIa mass & energy transfer off 


moftp-ftf at 


25 




No He reheating 




25 


100 


No SNII winds, no SNIa energy transfer 


NOSN N07rOOT 


25 


100 


^MTT winriQ nn S^TTn pnprov trnnQTPr pnrtlincT ii'ipc nnmorriinl nnnnrlnnnpc 


NOREION 


25 




No reionization 


NOZCOOL 


25 


100 


Cooling uses primordial abundances 




25 


100 






25 




RpHcriiTT rpi on 1 vciti on — f\ 


REIONZ12 


25 




Redshift reionization — 1 2 


cpAA/TPT y1 


25 






*iF*vT OPFlnTS 


25 




7ks = 1-75 


jr 1 rlKCjZj 


25 




Metallicity-dependent SF threshold 


cisJTor;ATTS*s 




100 


Gaussian SNIa delay distribution 


WriFMS 
vv i_jni>i o 


25 


100 


Wind mass loading and velocity determined by the local density 


WHYDRODEC 


25 




Wind particles temporarily hydrodynamically decoupled 


WML1V848 


25 


100 
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Figure 2. Metal distribution in the reference simulations. Shown are 5/i~^ Mpc (comoving) thick slices thi'ough the REFERENCEX025N5 12 simulation (top 
row) at z = 4 (left), and 2 = 2 (right) and through the REFERENCE_L 1 00n5 1 2 simulation (bottom row) at 2 = 2 (left), and 2 = (right). The colour coding 
shows (SPH-smoothed) metallicity, mass averaged along the line of sight; the colour scale is cut off below \og(Z / Zq) - -5 although metallicities extend to 
lower values. Red circles correspond to haloes identified by a friends-of-friends group finder, with radius proportional to the logarithm of the stellar mass of 
the halo, as indicated to the right of the panels. Metals are initially strongly clustered around haloes, but their volume filling factor increases as time progresses. 



© 0000 RAS, MNRAS 000, OOO-OOO 



6 R. P. C. Wiersma et al. 




Figure 3. Metallicity (top row), metal mass fraction (middle row), and Qz (bottom row) as a function of redshift for the various phases of baryons in the 
REFERENCEXIOOnS 12 simulation. Curves indicate simulation results, with different line styles indicating various baryonic phases. Left column: stars (solid 
black), SF gas (dashed red), NSF gas (dot-dashed blue); Middle column: diffuse IGM (black), cold halo gas (dashed red); Right column: WHIM (solid black 
line), ICM (dashed red). Data points indicate observatio ns and are colour-coded to indicate to which phase they should be com pared. The diamo nds al z = 2.5 
are estimates of stellar and star-forming metallicities by Pagel ( 2008), with further stellar metallicities indicated by a circle ( Gallazzi et al. 2008) and a square 
taalUdav et allEoOS). DLA measurements are from Prochaska et al. (2003, solid, left po inting triangles), IGM pixel optical depth measurements in QSO 
spectra are fr ormAguirre et alj j 2008t. open triangle pointed down) and lSchave et alj j2003l open triangle pointed up). The ICM measurements are from X-ray 
observations jSimionescu et alj2009l hatched region). 



predicted metallicities with observations and will mostly focus on 
metallicities in the rest of this work. 

Now we briefly walk through the data (which we have con- 
verted to our adopted solar abundances) to which we are compar- 
ing: 

• Stars For s tellar metallici t ies, w e compare to the global val- 
ues obtained b y iGallazzi et a l] ll2008l circle in the top-left panel) 
for 2 = and iHallidav et al.l ( l2008l square in top left panel) for 
z = 2. lHalliday et al. ( 2008) note that their metallicities are lower 
than t hose of Erb et al. ( 2006) and speculate that since lErb et al.] 
1 I2OO6I) observed oxygen, a non-solar [0/Fe] ratio might make their 
observations agree. We have also compared their point to [Fe/H] 
and indeed find a better match to the data. The 2 = 2.5 point is 
taken from lBouche et al] ( |2007|) which is compiled from their own 
previous measurements^ 

• Star-forming gas ISouche et alj ( |2007|) give a value for ISM 
and/or dust at z = 2.5 and, following iPagell we include it 
here as a filled red diamond in the top-left panel. 

• Diffuse IGM We have in cluded the carbon an d oxygen Lyman- 
a forest measurements from lSchave et al] (l2003l upwards pointed 
triangle in the top middle panel) and Aguirre et alj ( I2OO8I down- 
wards pointed triangle in the top middle panel), respectively. 

• Cold Halo Gas While the precise nature of DLAs is still some- 
what uncertain, we surmise that circum-halo and intrahalo cold gas 
most likely have higher cross-sections than the ISM and we thus 



assume that observations of DLA and sub-DLA systems trace cold 
halo gas. As such, we have p l otted t he compilation of DLA metal- 
licities from lProchaska et"al] ( l2003l leftwards pointed triangles in 
the top middle panel). 

• WHIM The WHIM has long stood as a phase that is thought 
to contain a large fraction of the baryons and metals, but is difficult 
to detect. As such, there is no direct measurement of the metallicity 
of the WHIM. 

• ICM Measurements of ICM metallicities are also trouble- 
some, although much less so than for the WHIM. The main dif- 
ficulty is measuring the metallicity out to large enough radii in or- 
der to get a good estimate of the mean metallicity. We have cho- 
sen to use the ou t ermos t measurements of cluster metallicities by 
ISimionescu et ah liooi hatched region in the top right panel), us- 
ing a box to indicate the range of values in their sample. 

The metallicities of our simulation follow the trends seen in 
the data. We reproduce the observed increase with cosmic time of 
the stellar and cold halo metallicities (although our stellar metal- 
licities fall slightly above the observed points). The high-redshift 
measurements for the star-forming gas are satisfied as well. The 
predicted metallicity of the 2 ~ 3 diffuse IGM agrees with the mea- 
surement for carbon (! Schave et al.l2003|)^ but is about 0.5 dex lower 
than the observed value for oxygen ( lAguirre et al. 2008), which 
should be closer to the overall metallicity. Our prediction for the 
intracluster metallicity is within the observed range. 
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Figure 4. Mean metallicity (left) and metal mass fraction (right) for stars (top), star-forming gas (middle), and non-star-forming gas (bottom) as a function of 
redshift for all OWLS simulations. Shown are the REFERENCE simulations (solid black), the rest of the l025n5 1 2 simulations (dotted red,) and the rest of the 
l100n5 12 simulations (dotted blue). The higher resolution (L025N5 12) simulations were stopped at z = 2 (with the exception of the REFERENCE model 
which was stopped just above z = 1), but the lower resolution simulations (L 1 00n5 1 2) were continued to z = 0. The mean stellar metallicity in units of solar 
rises from —0.5 dex at z = 6 to supersolar at z = 0, closely tracking the metallicity of the SF gas. For these two phases, resolution plays a minor role, with 
a 64 times higher mass resolution resulting in about 0.2 dex higher metallicities. The metallicity of NSF gas is, however, strongly resolution dependent for 
z > 2. In the higher resolution simul ation it rises from Lo g(Z) /Z(T) ^ —3 at z = 6 to to —1 at z = 0. Points with en'or bars indicate observations of stellar 
metallicities at z = and z = 2 bv lGallazzi et al.l j20O8h and lHallidav et alj )2008h , respectively. The stellar metal mass fraction increases towards z = 0, 
whereas the metal mass fraction in SF-gas decreases slowly to z i?i 2, then falls much more rapidly. Aside from the models without feedback (the group of 
vei7 low curves in the bottom panels; note that NOSN_l025n5 12 and NOSN_NOZCOOL_l025n5 12 fall on top of each other in the left panel and fall below 
the plotting area in the right panel), most of the different runs yield results that follow the REFERENCE model relatively closely, indicating that the predicted 
metallicities are reasonably robust to model variations. 



3.2 Overview 

We want to investigate which physical processes determine the cos- 
mic metal distribution in a given phase. We begin by discussing 
some of the overall trends that are present in nearly all of our 
models. Fig. |4] shows the evolution of the metallicity (left) and the 
metal mass fraction (right) for stars (top), SF gas (middle), and NSF 
gas (bottom). The black curves show the results for the reference 
model, with the curve that continues down to z = indicating the 
100 Mpc/h box. Blue and red curves indicate the results for the 100 
and 25 Mpc/h boxes for all other models. 

The metallicities of the SF gas and stars track each other very 



closely, increasing steadily from log{Z/ZQ) ~ —0.5 at 2 = 6 to 
supersolar values at jz = 0. While the metal mass fraction contained 
in stars increases with time, from ~ 10~^ at 2 = 6 to 0.2-0.9 at 
z — 0, the fraction of the metal mass contained in SF gas decreases 
with time, and does so more rapidly for 2 < 2. All models pre- 
dict that by 2 = SF gas accounts for only a small fraction of the 
metal mass. The rapid decrease in the metal mass fraction in SF 
gas for 2 < 2 reflects the decrease in the cosmic star formation 
rate jSchave et alj|2010l) and hence the global mass density of SF 
gas. While the metallicity of the NSF increases rapidly with time, 
from < 10"^ Z0 at 2: = 6 to 10"^ Zq at z ^ 0, the fraction of the 
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Figure 5. The left and middle panels show the fraction of metals in gas as a function of density and temperature, respectively, and the right panels shows the 
fraction of metals in stars as a function of stellar metallicity. The top and bottom rows show 2 = 2 and 0, respectively. Each histogram is normalised such that 
it sums to unity. The reference simulations are shown as solid black, while the rest of the simulations are shown in dotted red. The REFERENCE_l100n512 
model at z = 2 is shown in the top panel with the dashed blue line. SF gas is omitted from the central panel, since the 'temperature' of this multi-phase gas 
is set by the imposed P — p relation. The dotted vertical lines indicate the star formation threshold. While there are a few outliers, most models yield similar 
metal distributions, with the largest differences occuning at high p and small T. 



metals residing in this diffuse phase is relatively constant and, ex- 
cept for models without any feedback, always significant (for both 
reference simulations this fraction stays between 20% and 35% be- 
tween z = 6 and z — 1). Initially, the majority of the metals reside 
in SF gas. At z « 3 metals are approximately equally distributed 
over the three different phases (stars, SF and NSF gas), with stars 
becoming the dominant repository of metals at lower redshift. 

Comparing the two different box sizes (blue vs red curves), 
for which the mass resolution differs by nearly two orders of mag- 
nitude, we see that the results for stars and SF gas are surprisingly 
(given the enormous change in resolution) well converged. The 
higher resolution runs predict stellar metallicities and metal mass 
fractions that are about 0.2 dex higher. For SF gas the metallicity 
is also slightly higher (about a factor of two at 2; = 6 and less at 
lower redshifts), but the metal mass fraction is slightly lower, be- 
cause there is more gas in this phase in higher resolution runs. The 
metallicity of NSF gas is much more sensitive to the resolution and 
is unreliable in the lower resolution simulations above 2; ~ 2. The 
fact that the two different resolutions yield similar metal mass frac- 
tions for NSF gas at high redshift is a coincidence, caused by the 
fact that there are far fewer metals in the lower resolution runs. 

We recall the results of the more extensive numerical conver- 
genc e tests performed for t he reference model in appendices B and 
C of lWiersmaet^ ( I2OO9 I). For metallicity and metal mass frac- 
tion, except possibly for the ICM, our simulation boxes are suf- 
ficiently large. On the other hand, resolution proves to be much 
more of a challenge. The stellar metal mass fraction just barely 
converges for the l025n5 12 resolution, although the difference is 
small hy z — 2. The metal mass fraction in the cold-phase NSF gas 
is converged to within a factor of two. Obtaining converged results 
tends to be more challenging at higher redshifts. Metallicity is gen- 
erally better converged than the fraction of metals in a given phase. 



and is reasonably reliable for all phases at the resolution of the 
L 1 00n5 1 2 simulations, with the exception of the metallicity of the 
diffuse IGM. Higher resolution simulations generally yield higher 
metallicities, especially at higher redshift, but even in this phase, 
the simulations are converged by z = 2 (2: = 3) for l100n512 
(L025N512)runs. 

The distributions of metals at z = and 2 = 2 are investi- 
gated in more detail in Fig. [5] The metal mass weighted probability 
distribution function (PDF)|j shows several maxima in both density 
and temperature, with minima at p ~ 10'^ (p) and p ~ 10* at 2 = 2 
and 2 = respectively, which corresponds to nu ~ 10"'^ cm^'^ in 
both cases, and at T « lO'^ K. 

All models have a relatively well pronounced maximum at 
p ^ (p), which is due to the fact that there are very few metals 
in underdense gas and that gas that is much denser quickly accretes 
onto haloes. The outliers at the low-density end of the distribution 
result from turning off feedback (see Fig. |6}, while the outliers at 
the high-density end correspond to either strong feedback models 
(e.g. AGN) or models in which the star formation efficiency is in- 
creased (e.g. SFAMPLX3). 

The temperature minimum in the metal mass fraction at T ~ 
10^ K is caused by a number of factors. First, gas cooling is very 
efficient at this temperature. Second, the winds tend to shock met- 
als to temperatures higher than this value, causing a peak in the 
distribution at higher temperatures. Finally, the equilibrium tem- 
perature due to photo-heating is lower than this temperature, caus- 

We emphasise that the right most panel does not show the metallicity 
distribution function, as is commonly plotted. It rather shows how the metal 
mass is distributed among stars of various metallicities, which is more heav- 
ily weighted towards higher metallicities than the metallicity distribution 
function. 
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ing a second peak at T ~ 10* K. The temperature distribution of 
the metals shows less variation than the distribution with density. 
The two outliers again correspond to models without feedback (see 
Fig. |6l(. Such models have little shock-heated metals, resulting in 
much more low temperature metals. Without metal cooling, these 
metals pile up at 10* K, whereas with metal cooling the metals are 
found to very low temperatures. This effect is diminished at 2; = 2, 
although note that the NOSN_l025n5 12 simulation is not shown 
because it was stopped at z ~ 3 due to computational costs. 

The distribution of metals over stars with different metallici- 
ties (right panels) is similar in the different models, with the main 
deviation at high metallicity being due to the AGN simulation, 
which predicts more metal mass in lower metallicity stars. 

These figures once more illustrate the large dynamic range in 
density and temperature over which metals are distributed. These 
distributions are mostly converge d with respect to bo th box size 
and resolution (Figs. B3 and C3 o f lWiersma et al.l2009h . with some 
dependence of the temperature distribution on box size, and of the 
density distribution on resolution. Next, we will investigate in more 
detail how various physical processes affect the metal distribution. 

3.3 Impact of galactic winds and metal-line cooling 

In this section we discuss the effect of including galactic winds 
driven by SNe and metal-line cooling on the metal distribution. We 
show that outflows are essential for enriching the IGM and that 
metal-line cooling strongly affects the temperature distribution of 
metals. On the other hand, the hot non-star-forming gas can obtain 
rather high metallicities without feedback, indicating that the ICM 
is enriched at least partly by some process that is not directly related 
to SN feedback. 

Fig. |6] compares the REFERENCE simulation to a simulation 
without metal cooling (NOZCOOL), one without winds (NOSN; 
note that the blue lines in the lower left panel fall mainly under 
the green lines), and one with neither metal cooling nor winds 
(NOSN_nozcool). These simulations differ not only in the way 
metals are distributed, but also in their star formation histories since 
cooling and f eedback of course also affect the conversion of gas 
into stars (see lSchave et alj201(]|) . 

The top panels shows metal mass weighted PDFs at 2 = 0. 
The simulation without metal cooling contains fractionally more 
metals than the REFERENCE model at p ~ 10^ (p) — 10* (p) at the 
expense of metals in the ISM. Metal cooling thus allows metals to 
condense from haloes into star-forming gas. This clearly illustrates 
the importance of metal cooling for gas accretion onto star-forming 
regions. Surprisingly, the presence of metal cooling has little ef- 
fect on the temperature distribution where the cooling curve peaks 
(T ~ 10^ K). This is most likely because the cooling times here 
are already short. Metal cooling does reduce the fraction of met- 
als at T ~ lO'^'^ K. The largest difference is noticed in the low 
temperature regions (we remind the reader that the temperature 
distribution excludes gas with densities above our star-formation 
threshold). Metals can cool gas well below 10* K, while neglect- 
ing metal cooling abruptly cuts the distribution off, although there 
are still some metals with T < 10* K due to adiabatic cooling. Fi- 
nally, metal cooling slightly increases the fraction of metals locked 
up in high metallicity stars. 

In the absence of galactic winds, only a negligible fraction 
of the metals reach densities p^W (p) (top-left panel of Fig.|6ll. 
This underlines the importance of feedback in reproducing the met- 
als se en in the IGM. Note that this is in opposition to lOnedinl 
j 19981) who found that feedback played only a minor role in enrich- 



ing the IGM, but agrees well with lAguirre et al.l ( l200ll i%. Winds 
also serve to evacuate metals from the ISM as we can see that 
in the absence of feedback, metals tend to pile up just above the 
star-formation threshold. Without winds, the fraction of metals at 
(NOSN_NOZCOOL) or below (NOSN) 10* K is much higher than in 
any model that includes feedback from SNe (top-middle panel of 
Fig.lU, because the metals reside in gas with higher densities and 
metallicities and thus higher cooling rates. Winds strongly boost 
the fraction of metals in gas with 10^ K < T < lO'^ K, which 
suggests that wind shocks may also shape the metal distribution 
at these temperatures. As expected, galactic winds shift the metals 
locked up in stars to lower stellar metallicities. 

In the bottom panels of Fig.|6]we show the met allicities of var- 
ious p hases as a function of redshift. As discussed in lWiersma et al.l 
( l2009l) . in the REFERENCE model most of the metals are initially in 
the gas phase, while by the present day the majority of the met- 
als are locked up in stars. The effect of neglecting metal cooling 
is small for the metallicity of either stars or SF gas, but is obvious 
for the metallicity of the IGM (which is lower by up to an order of 
magnitude in the absence of metal cooling) and the WHIM (higher 
by ~ 0.2 dex without metal cooling). 

Without winds, more metals are in stars and significantly 
fewer in NSF gas, but the metallicities of stars, SF gas, dense NSF 
gas, and ICM are in fact not very different. The absence of feedback 
has a major effect on the metallicity of the IGM (low T and low p 
NSF gas) and the WHIM (hot, NSF gas), decreasing their metal- 
licities by more than 2 and 1-2 orders of magnitude, respectively. 
Without galactic winds, metals are not efficiently transported to the 
lower densities of the diffuse IGM or WHIM. A careful census of 
the IGM metallicity, and its evolution, is therefore a powerful probe 
of the properties of galactic winds through the ages. On the other 
hand, the metallicity of the ICM is hardly affected by the presence 
of winds, suggesting that dynamical processes such as ram pressure 
stripping dominate the enrichment of this phase. 

In summary, feedback from SNe is essential for enriching 
the gas outside of haloes and increase the fraction of the metals 
in warm-hot gas. Metal-line cooling, on the other hand, shifts the 
metal distribution towards lower temperatures and, provided winds 
are present, strongly increases the fraction of metals in the ISM. 
The metallicity of the present day ICM is ~ 10" * solar and is re- 
markably insensitive to the presence of feedback and metal cooling. 
This implies a minimum metallicity for the ICM, and an enrichment 
mechanism that is not related to galactic winds. 

3.4 Different implementations for galactic winds 

In the previous section we simply turned off galactic winds. Here 
we will investigate what happens if we vary the implementation of 
the outflows. We will show that the metallicity of the IGM is par- 
ticularly sensitive to the implementation of feedback, some mod- 
els are quite capable of enriching this low-density gas while others 
cannot. On the other hand, the metallicities of the stars and the star- 
forming gas are largely robust to differences in the wind model. The 
ICM and WHIM maintain a relatively constant metallicity with red- 
shift if the feedback recipe is varied. 

Winds in most of the simulations comprising the OWLS suite 
are characterised by two paramet ers: the mass loading factor, 77 , 
and the initial wind speed, jPalla Vecchia & Schavd 1200 8l : 
Eq. [3}. All models shown in Fig. |7] use the same value of the 
product 77 and hence they assume that the same fraction of 
core collapse SN energy is used to power winds. The wind speed 
sets an approximate maximum depth of the potential well out of 
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Figure 6. Distribution of metals for l100n512 simulations with different feedback and cooling prescriptions. Top panels: 2 = metal mass weighted PDFs 
of gas density, temperature, and stellar metallicity. Each PDF is normalised to unity. In the left-most panel, the star formation threshold (riu = 10~^ cm~^) is 
indicated by the dotted line. ISM gas was removed from the temperature PDF. Bottom panels : Metallicity as a function of redshift for the phases as defined in 
Fig. 1. The left panel shows stars (solid), SF gas (dashed) and NSF gas (dot-dashed). The middle panel shows cold IGM (NSF p < 10^ (p>, T < 10^ K; solid), 
and cold halo and circum-halo gas (NSF p > 10^ {p),T < 10^ K; dashed). Finally, the right panel shows the WHIM (NSF gas, 10^ K < T < lO''' K; solid) 
and the ICM (non-star-forming, T > 10^ K; dashed). Shown are the REFERENCE simulation (black), a simulation which ignored metal cooling (NOZCOOL; 
red), a simulation without galactic winds driven by SNe (NOSN; blue), and a model with neither winds nor metal cooHng (NOSN_NOZCOOL; green), as indicated 
in the top left panel. Note that the metallicity of the diffuse IGM for the simulations without feedback is below the plotted range in the bottom-middle panel. 
Data points indicate observations as in Fig. [3] where the points in the left panel coiTespond to stellar metallicities, the soHd triangles in the middle panel 
correspond to the cold halo gas while the open triangles show diffuse IGM measurements and the hatched region in the right panel shows an ICM metallicity 
measurement. Metal cooling shifts the metal distribution towards lower temperatures and increases the fraction of metals in the ISM. Galactic winds are crucial 
for the emichment of low-density gas. The metallicity of the ICM is insensitive to both metal cooling and galactic winds. 



which the wind can push baryons: star foimation will no longer be 
quenched significantly in a galaxy for which wind particles can- 
not escape. Note that since winds are quenched due to drag forces 
in the ISM, which depend on its pressure and hence on the mass 
of the galaxy, the relationship between the wind speed required 
for escape and the potential w ell of the halo is non-trivial (see 
iDalla Vecchia & Schavel I2OO8I) . This parametrisation of winds is 
still very simple with little physical motivation. Given the coarse 
resolution of cosmological simulations, wind models should be 
guided by observations. Unfortunately, relating the wind parame- 
ters to observations is problematic, because the wind mass load- 
ing and velocity will vary with radius (which is poorly constrained 
for measurements based on absorption lines), inclination, and gas 
phase. For a complete description of the method and further moti- 
vation see Dalla Vecchia & Schaye ( 2008). 

At high z, when most stars form in low-mass galaxies, even 
a low value of may enable the gas to escape from most galax- 
ies. For constant 7)u^, a low wind speed implies high mass load- 
ing and hence a strong suppression of early star formation. Metals 
are then very efficiently transported out of the galaxies, and end 
up predominantly in the cold NSF gas, since the wind velocities 
are low enough that strong shocks do not occur. However, as time 
progresses, a low wind speed model can no longer suppress star 
formation in the increasingly more massive galaxies. The metals 
that are produced can no longer escape from the haloes and remain 
in the SF gas. For higher values of metals can escape from typi- 



cal star-forming galaxies up to far lower z and end up in hotter gas 
because the winds get shocked to higher temperatures. 

The dependence of the z = 2 metal distribution on the wind 
speed and mass loading is shown in the top row of Fig.|7]for models 
that all use constant rj and Vw , identical riv^ , but for which -q varies 
by a factor of 8 (and thus by a factor of %/8). As expected, a 
higher Vw results in a larger fraction of metals residing in gas with 
low density and high temperature. Conversely, low-Huj models have 
a considerably larger fraction of their metals at high density, and 
low temperature (T ~ 10* K). The low mass-loading model shows 
a peak in the density distribution of the metal mass just above the 
star formation threshold (vertical dotted line) that is similar to the 
reference model, indicating that the location of this peak is not set 
by the wind velocity (although the width and the height certainly 
are). There is another small peak at p ~ 10^'^ (p) in the metal-mass 
weighted density PDF of the low mass-loading model, which may 
be due to the largest haloes which have such high ISM pressures 
that even these high-velocity winds are quenched. The metal mass 
distributions in temperature show how higher wind velocities lead 
to more metals in hotter gas, due to metals in winds getting shocked 
to higher temperatures. A large portion of metals that are at T ~ 
10^ K in the REFERENCE model are heated to > lO'' K for v^, = 
848 kms~^. The stellar metal mass PDF shifts to slightly lower 
metallicities for higher ■ 

While the distribution of the metals depends very strongly on 
the wind model, there is relatively little difference in the metallicity 
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Figure 7. As Fig.|6] but for l025n5 1 2 simulations with different values for the mass loading, r), and wind speed, v, but identical values for rjv^ and thus for 
the wind energy per unit stellar mass formed. The metal mass weighted PDFs (top row) are for z = 2. Compared are the REFERENCE simulation {r} = 2, 
= 600 kms~^; black), a simulation using a mass-loading of 1 (77 = 1, = 848 kms~^; red), and a simulation using a mass-loading of 8 (77 = 8, 
v-w = 300 kms~^; blue). The density and temperature distributions of metals are strongly dependent on 77, with higher values of rj shifting metals to higher 
gas densities and lower temperatures. Only the metallicity of the IGM is significantly affected, with higher rj yielding higher metallicities. 
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Figure 8. As Fig.|6] but for l025n5 12 simulations comparing the REFERENCE simulation (black) to model WDENS (red) in which wind speed is proportional 
to the local sound speed, while keeping the the wind energy per unit stellar mass formed identical do that in the REFERENCE model (Eq.|5). The metal mass 
weighted PDFs are for z = 2. The metallicities of both stars and SF gas are similar between these models, but the top-left panel shows that lowest density gas 
is significantly more enriched in the WDENS model. Nevertheless, the middle bottom panel shows that the metallicity of the cold IGM though is lower in the 
WDENS model, which is due to the fact that the metals are hotter in this model. 
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Figure 9. As Fig. |6] but for L025N512 simulations, comparing the REFERENCE simulation (black) to model WVCIRC (red), in which wind speed depends 
on the depth of the potential well (Eq. |6j, and to model WTHERMAL in which SN energy is injected thermally (blue). While WTHERMAL uses the same 
wind energy per unit stellar mass as the REFERENCE model, the energy increases with halo mass for model WVCIRC (and can exceed that available from 
star formation). The metal mass weighted PDFs are for 2 = 2. Model WVCIRC is able to evacuate metals from high-density regions, while at the same time, 
spreading the metals over a variety of temperatures. Model WTHERMAL differs less dramatically from the REFERENCE model, with more metals ending up in 
high-density regions. 



between models with different Vw , with the exception of the metal- 
licity of the IGM. In particular, the stellar and SF gas abundances 
are only marginally affected by the value of for the range of 
parameter values that we investigated. The metallicities of both the 
WHIM and the ICM remain strikingly constant with time at ~ — 1 
and ~ —0.5, respectively, all the way from z = 6 to z = 1 and, as 
shown in Fig.|6l even to z = 0. These values are clearly relatively 
robust with respect to the wind implementation. 

The IGM metallicity in the high-mass-loading (ri = 8) model 
is nearly an order of magnitude higher than in the low-mass-loading 
(ri — 1) model, despite the fact that a much smaller fraction of the 
metal mass resides at densities p < 10*^ (p) for rj = S. The increase 
in the metallicity of the IGM (which has T < 10'' K by definition) 
with decreasing must therefore be due to a decrease in the tem- 
perature of the enriched low-density gas. The rj — 8 model gives 
better agreement with the value inferred for the oxygen abundance 
bvlAg uirre et al. ( 2008), and is well within the uncertainties of the 
ISchave et alj ( |2003|) value for carbon. Because most of the metals 
in the IGM are ejecte d by low-mass galaxies (see lBooth et al.l2010l : 
Iwiersma et al.ll201oh from which even low-velocity winds can re- 
move gas, a large mass-loading enriches the low-density IGM much 
more efficiently. 

In model WDENS the fraction of SN energy that powers the 
wind, esN oc rj (Eq. [3]l, is still kept fixed to the value used in 
the reference model, but scales with the local sound speed, Cs- 
Since star-forming gas is on an imposed P ~ p relation (Eq. [T](, 
Ca oc P/ p oc pT<=« ~^ oc p^^^ for our assumed value of 7cfi — 4/3. 
The assumed wind parameters therefore depend on the density: 



?^H:thr 



1/6 



and 



(4) 



-2/6 



(5) 



where nH;thrcs is the star formation threshold density. We chose 
vo = 600 km s^^ and -qo — 2, so that the wind parameters are the 
same as in the REFERENCE model at the star formation threshold. 
Compared with the REFERENCE model, the winds in WDENS will 
remain efficient in regulating star formation for denser gas at the 
bottom of deeper potential wells, where the winds will have higher 

Figure [8] confirms that model WDENS is able to reduce the 
build up of metals in high-density gas and the fraction of metals 
locked up in high-metallicity stars. In general, WDENS continues 
the trend with increasing shown in Fig. |7] The high-density 
tail is truncated at p/{p) ~ 10^. The higher wind speeds result 
in a larger fraction of metals at T > 10^ K due to winds getting 
shocked. With fewer metals in high-density SF gas, stellar metallic- 
ities are on average « 0. 1 dex lower than in the REFERENCE model. 
The ICM is enriched to higher levels at nearly solar, with little evo- 
lution. We note, however, that at these high redshifts the 'ICM' 
consists largely of hot, low-density gas shocked by galactic winds. 
As we will show below (see Fig. |13t , for 2 = model WDENS 
predicts ICM metallicities that are only slightly greater than for the 
reference model. The enrichment of the cold IGM is notably less 
advanced than in the REFERENCE model, with a metallicity nearly 
1 dex lower, in spite of the larger fraction of metals at low densi- 
ties. This is again because the faster winds shock-heat more of the 
escaping metals to temperatures for which the cooling time is long. 

In model WDENS the wind speed scales with the local sound 
speed to obtain efficient feedback in low-mass galaxies at early 
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timelfl as well as in more massive haloes at later times. An alterna- 
tive method to obtain a similar behaviour makes the wind speed de- 
pend dire ctly on the depth of the potenti al well, as in the model dis- 
cussed bv lOppenheimer & Pavel fzOOSi) . which is meant to mimic 
a case where the wind is driven by radiation pressure. Such a wind 
is implemented in our WVCIRC model as follows: an on-the-fly 
friends-of-friends group finding algorithm is used to identify haloes 
and their circular velocities, Wciic ~ \/ GMjR^i^ are computed. 
The wind speed and mass loading are then given by 



(3 + r?o ) 



-«cir 



(6) 



(7) 



where Ucrit = 150 kms"'^ and rjo — 2 dOppenheimer & Davd 
[2008). In this model the effective value of egjv is no longer held 
constant, as the momentum carried by the wind is kept fixed in- 
stead. For galaxies in groups and clusters the wind energy exceeds 
that available from SNe. This is not necessarily a problem, as the 
winds may be d riven by radiation pre s sure. However, for the param- 
eter values of Oppenhei mer & Pavel ( I2008 1. the momentum that is 
injected also exceeds that available in the form of radiation (Haas 
et al., in preparation). These caveats should be kept in mind. 

Figure |9| shows that model WVCIRC transports significantly 
more metals to low densities, strongly reducing the fraction of met- 
als located up in high-metallicity stars, similar to WDENS (Fig.jSj. 
However the temperature distribution of the metals in this model is 
quite different from that of WDENS with a much smaller fraction of 
metals at ~ 10® K and a much larger fraction at ~ 10^ K. WVCIRC 
predicts a metallicity of the diffuse IGM that is up to an order of 
magnitude higher than for the REFERENCE model. This is a con- 
sequence of the high mass loading factors used in small galaxies, 
which results in much more efficient transport from galaxies to the 
IGM, combined with the low wind velocities, which ensure that 
the winds do not shock the gas to high temperatures. The higher 
metallicity of the cold IGM ag rees better with the ob servations, as 
was already shown by O ppenheimer & Pavel ^20061) for a similar 
model. Stellar metallicities, which are a factor of « 2 lower by 
~ 2, also agree better with observations. 

In Fig. [9] we al so show resul t s from another model, WTHER- 
MAL, discussed in ISchave et al.l feOld) and Palla Vecchia and 
Schaye (in prep.), in which the energy that is injected is again iden- 
tical to that in the REFERENCE model, but in which this SN energy 
is deposited locally as thermal rather than kinetic energy. As star 
formation occurs in high density regions that cool efficiently, such 
a feedback mechanism may be very i nefficient i f the multi-phase 
nature of the ISM is not resolved (e.g. lKatj|l992l) . Our implemen- 
tation of thermal feedback avoids much of the overcooling prob- 
lem by giving the energy to only a few particles so that they can 
be heated to s ufficiently high temp eratures that the cooling time is 
long (see also iTheuns et al.ll200a) . Nevertheless, model WTHER- 
MAL predicts a greater fraction of metals to reside in dense gas and 
high-metallicity stars, indicating that the winds in this model are 
less efficient than in the REFERENCE model. The metallicity of the 
cold IGM is, however, higher, because of the lower temperature of 
the enriched gas, which is no longer shock-heated by high-velocity 
winds. Consequently, the metallicity of the WHIM is slightly lower 
for model WTHERMAL. 



^ Note that lowering vq, and thus increasing rjo, would further boost the 
efficiency of the feedback in low-mass galaxies. 



In summary, the details of the wind implementation have a 
large impact on the predicted metal census. In particular, the frac- 
tion of metals in a given phase varies strongly with the initial wind 
velocity and mass loading, even if the energy in the wind is kept 
fixed. Compared with the large range in metallicity spanned by the 
various phases, the effects of variations in the wind implementa- 
tion on the metallicities are mostly minor. The metallicity of the 
cold IGM is most sensitive to wind implementation. While mod- 
els with high wind velocities, and thus low mass loading factors, 
shift metals to lower densities, they tend to predict lower metal- 
licities for the cold IGM, as the metals are shocked to high tem- 
peratures. Indeed, the high observed metallicities for this phase are 
better matched by models with low wind velocities. High wind ve- 
locities are, however, required for at least the more massive galaxies 
to reduce the stellar metallicities to the observed values. If we allow 
the wind parameters and energy to vary with the mass of the halo, 
with high mass loading factors in low-mass haloes and high wind 
velocities in high-mass haloes, then we can simultaneously boost 
the metallicity of the cold IGM and reduce that of the stars. This is 
the case for model WVCIRC( Oppenheimer & Pa vdl2008l) . although 
we noted that this model uses both more momentum and energy 
than is available from star formation. Finally, we note that it would 
be worthwhile to analyse the gas abundances in the simulations in 
a way that better mimics the observational determinations. 



3.5 Active Galactic Nuclei 

Feedback from accreting black holes in Active Galactic Nuclei 
(AGN) has for example been invoked to reproduce the observed 
black hole scaling relations, the break in the galaxy luminosity 
function, and the thermodynamic properties of groups and clus- 
ters of galaxies (e.g. Silk & Rees 1998; Kauf fmann & Haehneltl 

3d; 



l2000lJPi Matteo et al ,.,2005;. Bower et al..2006.:^c'roton et alj2006 
Booth & Schaye 200^ l2010| - ElcCarthv et al.ll2010l) . 

Our models for black hole formation and the associ- 
ated AGN feedbac k whi ch are modified versions of the model 
of [Sp ringel et aL (|2005|) are fully describe d and tested in 
iBooth & Schaveil2009h (see ISchave et al.ll20I(]| for a comprehen- 
sive summary) and we therefore only touch on the details here. We 
seed low-mass haloes with central black holes, which are then al- 
lowed to grow by mergers and gas accretion. A fixed fraction of 
1.5% of the accreted rest mass energy is injected into the ambi- 
ent gas by h eating it to a very high temperature. This model re- 
produces the lMagorrian et aU ( 1 19981) relation, the normalisation of 



which is determined by the assumed efficiency factor of the feed- 
back, as well as the other o bserved black hole scaling relations 
teooth & Schavel2009l , l2010l) and bo th the optical and X-ra y prop- 
erties of observed groups of galaxies l lMcCarthvet"aLll2010l) . 

Fig. [ini shows that AGN feedback is able to remove metals 
from high-density regions and transport them outside haloes (i.e., 
p < 200 (p)), while heating them to high temperatures. The metal- 
licity of the stars that contain most of the metals is reduced by 
nearly an order of magnitude. The metallicity evolution of the other 
phases is similar to the REFERENCE model, although the metallic- 
ity of the ICM is about a factor of two higher. Since there are less 
metals in the star-forming regions, it is not surprising that model 
AGN has lower stellar metallicities. However the WHIM's metal- 
licity is also slightly reduced, notwithstanding the larger fraction 
of metals in the WHIM (as seen in the temperature PDF). AGN 
feedback (as implemented here) therefore increases the fraction of 
baryons in the WHIM. 
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Figure 10. As Fig. |6] but for L100N512 simulations comparing the REFERENCE simulation (black) to model AGN (red), which includes feedback from 
accreting black holes. The metal mass weighted PDFs are for z = 0. AGN feedback is able to evacuate metals from the high-density regions and decreases 
the metallicity of stars at 2 = by 0.5 dex. The metallicities of most of the other phases remain largely unchanged with the exception of the hottest phase 
whose metallicity increases by a about a factor of two. 
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Figure 11. As Fig.|6] but for l10()n5 12 simulations comparing the REFERENCE simulation (black), with two models with a top-heavy IMF at high-pressure; 
one in which the extra wind energy is used to increase the wind velocity (ti^, = 1618 km s~^, DBLIMFC0NTSFV161 8,red), and one in which the extra wind 
energy is used to increase the mass loading {rj = 14, DBLIMFCONTSFML14, blue). The metal mass weighted PDFs are for 2 = 0. In the top-left panel, the 
black and red dotted lines indicate the star formation threshold and the density above which the IMF becomes top-heavy, respectively. A top-heavy IMF in 
starbursts only makes a difference if the extra SN energy is used to increase the wind velocity, in which case it shifts metals to lower densities and higher 
densities and reduces the stellar metallicity. 
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3.6 A top-heavy IMF in starbursts 

We now consider the effect of changing the IMF on the cosmic 
metal distribution. We will show that employing a top-heavy IMF 
in high-pressure regions affects mainly the stellar metallicity, leav- 
ing the metallicities of the other phases mostly unchanged. This 
is largely a result of the increased energy released from the SNe 
which are more prevalent for a top-heavy IMF. 

There is tentative observational evidence that the IMF in some 
environments is 'top-heavy' as compared to the solar neighbour- 
hood, such as in intense starb ursts ( McCradv et al. 2003, e.g.), or 
the inner parsec of the Galaxy dManess et al.l2007l . e.g.). Models of 
galaxy formation hint that such a top-heavy IMF may be required 
to ex plain galaxy counts in the sub mm and far IR jLacev et al.l 
I2OO8I e.g.). fPadoan et all ( Il997l) and lLarsorJ J2005l) discuss theo- 
retical moti vations for environm ental dependence and evolution of 
the IMF, see lLacev et alj ( l2008h for more discussion. 

We have implemented an environment-dependent IMF in our 
simulations, labelled DBLIMF in Table [T] in which the IMF is 
changed from the Chabrier fit to become top-heavy, $(M) oc M^^ 
for all M (compared with $(A/) oc Kr'^'^ for M > IMq for 
Chabrier), above a density threshold of nn ~ 30 particles cm~^, 
which corresponds to a pressure of P/k — 2.0 x 10® cm"^ K. These 
values were chosen because they result in ~ 10% of the stellar mass 
forming with a top-heavy IMF. Because the star formation rate den- 
sity increases with pressure, the IMF becomes top-heavy in regions 
with high star formation densities, i.e. 'starbursts'. For the sake of 
brevity, we restrict ourselves to the case in which the star forma- 
tion law is continuous across the 'starburst' pressure threshold (see 
[Schave et alj20ld for a full discussion of this point; see section lTTl 
for the results of models in which this is not the case). Note that 
these models are likely less well converged with resp ect to the nu- 
merical resolution than the REFERENCE model (see ISchave et al.l 
[2OIG). 

A top-heavy IMF increases the number of core-collapse SNe 
per unit stellar mass. We account for this by increasing the energy 
per unit stellar mass that is injected into the winds. We explore two 
extremes, one uses a higher mass-loading but the standard wind 
speed, (?7 = 14, Hn, — 600 kms^^), the other uses the standard 
mass loading but higher wind speed, {j] = 2,Vw = 1618 kms~^); 
note that this change is only enacted above the density threshold 
for top-heavy star formation. 

In Fig.[TT]we compare DBLIMFCONTSFV1618 (77 = 2,v^ = 
1648 kms"^) and DBLIMFCONTSFML14 (77 = 14, = 
600 km s"^) with the REFERENCE model. While a top-heavy IMF 
increases the amount of metals produced per unit stellar mass in 
high-density regions due to the enhanced fraction of high-mass 
(high yield, short lifetime) stars, the increased feedback serves to 
move a lot of these metals to high temperatures and low densities, 
resulting in a slightly lower stellar metallicity, but a higher ICM 
metallicity for DBLIMFCONTSFVI6I8. Channeling the extra feed- 
back energy into a higher mass loading factor yields metallicities 
and metal mass distributions that are nearly identical to those of the 
REFERENCE model. This is because feedback in high-pressure gas 
is inefficient for such a low wind velocity. 

The use of a top-heavy IMF in 'starbursts', as implemented 
in the DBLIMF models, has surprisingly little effect on the relative 
abundance of Q-enriched elements produced by these massive stars 
(not plotted). While [O/Fe] averaged over entire simulation volume 
is higher for a top-heavy IMF, the diffuse IGM and cold halo gas 
have a lower a enhancement if the IMF is top heavy, because the 
a elements are deposited in hotter gas. This confirms that the dif- 



ferences between the REFERENCE model and the top-heavy IMF in 
starbursts models are driven by the changes in the feedback rather 
than by the increased metal production. 



3.7 Other models 

Although some processes can dramatically affect the cosmic metal 
distribution, we will see in Figs. ll2lfT5] that there are a number of 
models that show little difference. It is remarkable that many dras- 
tic changes in physical parameters and processes have little affect 
on the metallicity and metal mass distribution of most, or even all, 
phases. We briefly comment on each of the models that exhibit this 
behaviour and refer to Schave et al. (2010) for more detailed de- 
scriptions of the models. 

NOREION, REIONZ06, REIONZ12 - We have also varied 
the reionization history of hydrogen in our simulations. Whereas 
the REFERENCE model assume H reionization at Zr — 9, these 
models have no reionization, Zr — 6 and 12, respectively. We find 
that so long as the heating from reionization happens, the simula- 
tions have a very 'short' memory of their thermal history. If there 
is no heating from reionization, then a slightly larger fraction of the 
metals lie in cold gas, although the metallicity of cold, dense gas is 
actually a little lower. 

Aside from directly heating the mass, the reionization history 
determines the lowest halo mass in which galaxi es are able to f orm. 
This may be espe cially important since Wiersm a et al. and 
iBooth et all ^2010 ) showed that the diffuse IGM is typically en- 
riched by low-mass haloes at high redshift. Between z — 9 and 
2 = 0, the virial mass of corresponding to 10* K increases from 
2 X 10* M0 to « 3 X lO'-* M0. Even the latter corresponds to 
only ~ 3 X 10^ dark matter particles in our 25h~^ Mpc box. It 
is therefore likely that we have strongly underestimated the impor- 
tance of the reionisation history for the metal content of the IGM. 

SFAMPLx3, SFSLOPElp75, SFTHRESZ - The normalisa- 
tion, power-law index, and threshold gas density of the Kennicutt- 
Schmidt star formation law (Eq.O are varied. The first model as- 
sumes a normalisation that is higher by a factor of 3; the second 
model has a slope of 1.75 (as opposed to 1.4); in the last model the 
threshold for star formation decreases with increasing metallicity. 
For models with more efficient star formation, the fraction of the 
metals residing in the ISM is lower. This is because feedback allows 
star formation to self-regulate, such that the rate of energy injection 
is independent of the star formation law. Hence, a shorter gas con- 
sumption time scale is balanced by the ejection of gas (and metals) 
from galaxies, su ch that the star formation rate remains the same 
dSchave et ai]|201Q) . As the ISM has a higher metallicity than the 
gas at larger distances, this transfer of metals increases the metal- 
licity of cold, halo gas, while leaving the metallicity of the ISM 
unchanged. 

NOAGBJVOSNIa - The contribution from A GB and type 
la SN e to stellar mass loss and feedback is ignored. 1 Schave et al.l 
(l2010h showed that this change reduces the z = star formation 
rate by about 40%. We find that the metallicities drop by a simi- 
lar factor, but that the metal mass fractions remain essentially un- 
changed. The abundance ratios (not shown) evolve very little and 
reflect solely the SNII yields ([0/Fe]^ 0.7 for all redshifts and 
phases). 

EOSlpO, EOSlp67 - The slope, 7off, of the p cc relation 
imposed on star-forming gas, is changed from the value jcS = 4/3 
assumed in the REFERENCE model to 1 and 5/3, res pectively. As 
was the case for the cosmic star formation history jSchave et al.l 
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bold) , this change has no significant effect on the metal distribu- 
tion. 

IMFSALP - Assumes a Salpeter IMF, as opposed to the 
Chabrier IMF. The change in IMF results in a small (less than a fac- 
tor of two) decrease in the metallicities and [0/Fe] ratios (not plot- 
ted) of all phases, which reflects the lower rate of metal production 
per unit stellar mass and the lower star formation ra tes that result 
from the associated reduction in the cooling rates (see lSchave et al.l 
[20l3). 

WHYDRODEC - Wind particles are temporarily 'decoupled' 
from the hydr odynamics (i.e. they feel only gravity) as in the widely 
used model o f ISpringel & Hemguistl ( l2003h . This modification has 
a significant effect, decreasing the metallicities of stars and the cold 
IGM by about a factor of two. Moreover, it roughly doubles the 
fraction of metals in, but not the metallicity of, the WHIM at the 
expense of the fraction locked up in stars. We have relegated this 
simulation to this overview due to space constraints, but it is an 
important variation. 

MILL - This mo del use the cosm ological parameters of the 
Millennium simulation jSpringellZOOSl) . as opposed to the WMAP3 
values. The most significant difference is the larger value of erg in 
the MILL model (0.9 vs. 0.74), which implies that structure forms 
slightly earlier. The wind mass loading is doubled (r) — 4:,Vw = 
600 kms^^) in order to r etain the match to the peak of the cos- 
mic star formation history jSchave et alj2O10l) . Compared with the 
REFERENCE model, all of the curves shift to higher redshifts, but 
as Figs. 11211151 show, even such a large change in cosmology has 
only a minor effect on the predictions for z — and 2. 

DBLIMFML14, DBLIMFV1618 - These models assume a 
top-heavy IMF for stars formed in high-pressure regions. They are 
identical to the models we discussed in section [X6l except that the 
star formation law is discontinuous across the 'starburst' thresh- 
old. Since no discontinuities in star formation density (as a func- 
tion of gas density) have been observed, this model decreases the 
efficiency of star formation in high-pressure regions in order to pro- 
duce a continuous UV flux as a function of pressure. The fact that 
these models show little difference to their DBLIMFCONTSF coun- 
terparts is consistent with our finding that variations in the star for- 
mation law are unimportant. 

NOHeHEAT - As discussed in section [2] we include some 
heating around z ~ 3 to account for the rise in the IGM tempera- 
ture around the time of helium reionization. Ignoring this extra heat 
input does not yield any significan t differences, as wa s also true for 
the cosmic star formation history dSchave et al]|20IOl) . 

SNIaGAUS - Given the large uncertainties in the normali- 
sation and shape of the cosmic SNIa rate, we have run a model 
where it is based on a Gaussian delay function (rather than the e- 
folding delay function used in the REFERENCE model). This results 
in roughly the same number of SNIa, but they happen much later in 
the history of the universe. The evolution of the [0/Fe] ratio shows 
the expected differences (the dip to lower values occurring at lower 
redshift), but the overall distribution of metals and metallicity are 
nearly identical to the results for the REFERENCE model. 
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Figure 12. Fraction of metal mass in various phases for a selection of 
L 1 00n5 1 2 models of the OWLS suite at ^ = 0. The orange, green, magenta, 
blue and red bars coirespond to the metal mass fraction in star-fomiing gas, 
ICM, WHIM, cold halo gas, and diffuse IGM, respectively. The remain- 
der of metals is locked up in stars. Metal fractions depend largely on the 
strength of feedback - models without feedback have very high stellar metal 
fractions, whereas models with strong feedback predict that most metals re- 
side in the WHIM. 
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3.8 Summary 

Figures [72l - 1 1 51 summarise our different models, comparing metal 
fractions and metallicities at redshifts z = for the lOO/i^^ Mpc 
and at 2: = 2 for the 25h^^ Mpc simulations. Note that many mod- 
els were only run in one of the two boxes (see Table [Hi. We give 
observational estimates of metallicities to compare to, but we cau- 



Figure 13. Metallicities of various phases for a selection of l100n512 
models of the OWLS suite at 2 = 0. The values for the REFERENCE model 
are in dicated by the horizo ntal lin es. Observations shown are as follows: 
stars - iGallazzi et"ai ? llOOS); ICM - ISimionescu et alJi2009l) : the rest - esti- 
mates compiled by Paget 1 200§). The metallicities are surprisingly similar 
across most models, but models without feedback are clearly different. 
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efficient feedback, which result in much higher metal fractions for 
the WHIM, at the expense of that of the stellar component. While 
the stellar metallicities are also lower in these models, the WHIM 
has nearly the same metallicity, indicating that the change in its 
metal mass fraction is mostly due to a change in its baryonic con- 
tent. 

The data at z = 2 are more uncertain and more incomplete 
(Figs. [141 andllSt. The models that stood out at z = are also very 
different at this redshift. Models with a top-heavy IMF and high 
wind velocity, the simulations with hydrodynamically decoupled 
winds, and model WDENS are generally similar to AGN and WV- 
CIRC, although the latter stands out in terms of the high metallicity 
and metal fraction that it predicts for the cold IGM. 
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Figure 14. As Fig.[T2] but for z = 2. 
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Figure 15. As Fig.|J3l but for z = 2. Observations shown are as follows: 



stars - lHallida_ 



vet al. 12008), 



cold halo gas - Prochaska et alj and 
diffuse IGM - lAguirre et al., (2008i). The diffuse IGM shows the most vari- 
ation among the diiferent models. 



tion that these are only meant to give an idea as to what values are 
reasonable, as the uncertainties are large. 

At redshift z — 0, the models that stand out most in 
terms of the metal mass fractions are AGN, WVCIRC, NOSN, and 
NOSN_NOZCOOL (Fig.[T2l(. The NOSN models (no energy feedback 
from SNe) predict significantly higher metal fractions in stars, and 
also higher stellar and ISM metallicities. Fig. 1 131 shows that these 
models also stand out in terms of the low metallicity of the WHIM 
and the diffuse IGM. Models AGN and WVCIRC include the most 



4 CONCLUSIONS 

We have investigated the distribution of metals predicted by simu- 
lations ta ken from the Overw helmingly Large Simulations project 
(OWLS: ISchave et al.ll2010 [). We considered different recipes and 
physical parameters in order to determine what the factors are that 
shape the distribution of metals as a function of the gas density, 
temperature, and metallicity. This builds on our previous work 
where we introduced our method and demonstrated th e effect of 
simulation box size and resolution dWiersma et alj200 9'). 

Our results can be divided between statements about the ef- 
fects of physical processes on the metal mass distribution and state- 
ments about the metal mass distribution itself. The former can be 
summarised as follows: 

• Metal-line cooling increases the fraction of the metals that re- 
side in stars and the ISM, at the expense of the fraction in highly 
overdense warm-hot gas. Metal-line cooling increases the metallic- 
ities of all cold phases. In particular, it increases the metallicities of 
the diffuse IGM and cold halo gas by about a factor of two. 

• While other mechanisms can enrich the warm-cold IGM (e.g. 
ram pressure and tidal stripping), without SN feedback, its metal- 
licity is far smaller and is strongly ruled out by observations of 
QSO absorption lines. 

• Processes other than outflows driven by feedback from star 
formation and AGN contribute significantly to, or even dominate, 
the enrichment of the ICM. 

• Even for a constant energy per unit stellar mass formed (i.e. 
a fixed feedback efficiency), the freedom given by the parameters 
of sub-grid implementations of galactic winds is very large. At a 
fixed efficiency, higher wind velocities (and thus correspondingly 
low mass loading factors) increase the fraction of the metals that 
reside in low-density (p < 10^ (p)) gas. However, since higher 
wind velocities also shift metals to higher temperatures, the metal- 
licity of the diffuse, warm-cold IGM is in fact higher for lower wind 
velocities. 

• The metallicity of the warm-cold diffuse IGM is most sensi- 
tive to the implementation of galactic winds. 

• By varying the parameters of the wind model with local prop- 
erties, one can significantly affect the metal distribution in the IGM. 
For example, a model in which low-mass galaxies drive highly 
mass-loaded, slow winds while high- mass galaxies drive fast winds 
with low mass loading factors (e.g. lOppenheimer & Davdl2008l) 
can efficiently enrich the diffuse warm-cold IGM while still limit- 
ing the build up of metals in the stellar components of more massive 
galaxies. 

• Feedback from AGN strongly reduces the fractions of metals 
that reside in the ISM and stars as well as the metallicities of these 
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two components. Fast galactic winds in massive galaxies, e.g. re- 
sulting from a top-heavy IMF at high gas pressures, work in the 
same direction but have a smaller effect than AGN. Efficient feed- 
back in high-mass galaxies only has a relatively small effect on the 
metallicities of the diffuse gas components. 

• Changes in the star formation law, the structure of the ISM, 
the cosmology, and the IMF all play a role, but are less important 
than the inclusion and implementation of metal-line cooling and, 
most importantly, outflows driven by feedback from star formation 
and accreting supermassive black holes. We also find a small effect 
for changes in the reionisation history, but for the case of the IGM 
this may be due to our Umited resolution. 

We summarise the metal mass distribution as follows: 

• Stars and the WHIM are the dominant depositories of metals. 
In all models they contain together at least 78 % of the metals at 
a = and at least 53 % at z = 2. At high redshift the ISM also 
contains a large fraction of the metals: at least 1 1 % in all simula- 
tions at z = 2. In all our models, and at both z = and z = 2, 
the remaining components (diffuse cold- warm IGM, cold halo gas, 
and ICM) together contain less than a quarter of the metals. 

• The mean metallicities of the WHIM and ICM are nearly con- 
stant in time. In most models they increase slightly with time to 
~ 10~^ Zq at z = 0. In contrast, the metallicities of the cold halo 
gas and the diffuse IGM, increase by more than an order of magni- 
tude from z = 5 to 0. 

• All our models predict, both for « = and 2, mean metallici- 
ties of Z ~ for stars and ISM, and 10"^ Z,j < Z < Zq for the 
ICM and for cold halo gas. Except for the models without SN feed- 
back, all simulations predict, both for z = and 2, mean metallici- 
ties of Z ~ 10"^ Zq for the WHIM and 10"^ < Z < 10"^ Zq 
for the diffuse cold-warm IGM. 

The inclusion and implementation of outflows driven by feed- 
back from star formation and AGN is most important for predic- 
tions of the distribution of metals. As cosmological simulations 
will need to continue to use sub-grid implementations for these 
processes for the foreseeable future, this implies that predicting 
the distribution of metals will remain difficult for ab initio mod- 
els. However, as we have shown, provided that some feedback is 
included, the mean metallicities of most different components can 
be robustly predicted to order of magnitude. 

The sensitivity of some of the results to feedback can also be 
regarded as fortunate. It gives us a useful tool to study the effects of 
outflows, which are an essential but very poorly understood ingre- 
dient of models of galaxy formation and evolution. The metallic- 
ity of the diffuse warm-cold IGM, which can be constrained using 
QSO absorption lines, is particularly sensitive to variations in the 
models and is therefore a very promising probe of the physics of 
galactic outflows. 
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